Electron-phonon interaction in Strongly Correlated Systems 
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The Hubbard-Holstein model is a simple model including both electron-phonon interaction and 
electron-electron correlations. We review a body of theoretical work investigating the effects of 
strong correlations on the electron-phonon interaction. We focus on the regime, relevant to high-T c 
superconductors, in which the electron correlations are dominant. We find that electron-phonon 
interaction can still have important signatures, even if many anomalies appear, and the overall 
effect is far from conventional. In particular in the paramagnetic phase the effects of phonons are 
much reduced in the low-energy properties, while the high-energy physics can be strongly affected by 
phonons. Moreover, the electron-phonon interaction can still give rise to important effects, like phase 
separation and charge-ordering, and it assumes a predominance of forward scattering even if the 
bare interaction is assumed to be local (momentum independent). Antiferromagnetic correlations 
reduce the screening effects due to electron-electron interactions and revive the electron-phonon 
effects. 

PACS numbers: 71.27.+a, 71.10.Fd, 71.30.+h, 71.38.-k 
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INTRODUCTION 



A wealth of materials, including the most challeng- 
ing systems (cuprates, manganites, fullerenes, ...), 
present clear signatures of both electron-electron (e-e) 
and electron-phonon (e-ph) interactions, leading to a 



competition -or interplay- which can give rise to differ- 
ent physics according to the value of relevant control pa- 
rameters and of the chemical and electronic properties 
of the materials. In this topical review we will focus 
on high-temperature superconductors, with the copper- 
oxide compounds (cuprates) in a prominent role, and an 
attention to the alkali-doped fullerides. 

In the case of the cuprates, which are arguably the 
most accurately studied materials in the last twenty-five 
years, the signatures of electron-phonon interactions are 
nowadays clear [lj, eventhough the overall scenario is far 
from ordinary: Electron-phonon fingerprints are evident 
in some properties, while they are weak or absent in 
other observables. Specifically, clear polaronic features 
are observed in optical conductivity [2J as well as in angle- 
resolved photoemission experiments (ARPES)Q in very 
lightly doped compounds. A substantial e-ph coupling 
can also be inferred by the Fano line shapes of phonons 
in Raman spectra and by the rather large frequency 
shift and linewidth broadening of some phonons at T c . 
Phonons are also good candidates to account for the 
famous "kink" in the electronic dispersions observed in 
ARPES experiments On the other hand phonons, 

which typically affect resistivity in standard metals, do 
not appear in transport experiments on cuprates. For 
instance the resistivity around optimal doping is ubiqui- 
tously linear in temperature (even in systems with rela- 
tively low critical temperature) and no high-temperature 
saturation seems to be present up to the highest achieved 
temperatures. This dichotomic behavior of cuprates, 
which display clear phononic features in some cases and 
none in others is one of the puzzling and challenging is- 
sues raised in these materials. Yet, it is almost univer- 
sally recognized that the main player in the cuprate game 
is the e-e correlation, which drives the parent compound 
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insulating and it is probably responsible of the main fea- 
tures of the overall phase diagram, including d-wave su- 
perconductivity and the existence of a pseudogap above 
the superconducting critical temperature. Therefore, it 
is not surprising that the signatures of e-ph interaction 
in the cuprates can hardly be understood in terms of 
the standard theory of e-ph interactions in weakly corre- 
lated metals, and a new theoretical framework including 
e-e correlations is needed. We will argue here that this 
change of perspective can indeed reconcile the different 
relevance of phonons in the various observables in corre- 
lated systems. 

On the other hand the superconducting members of the 
fulleride family, of composition A3C60 with A an alkali- 
metal atom, are often considered as standard phononic 
superconductors, in which the coupling between electrons 
and the local vibrations associated to the distortion of the 
carbon bukyballs is the driving force of superconducting 
pairing. The conventional nature of these compounds is 
challenged by recent investigations in expanded fullerides 
revealing several physical properties associated to e-e cor- 
relations. Indeed the CS3C60 solid with A15 structure is 
an antiferromagnetic Mott insulator at ambient pressure 
which becomes superconducting only under applied pres- 
sure, with T c reaching 38K0. The phase diagram as a 
function of pressure closely resembles that of cuprates as 
a function of doping, suggesting a central role of correla- 
tions. Consequently, e-e interactions are expected to be 
important in other members of the fulleride family. In- 
deed it has been shown that, thanks to the orbital degen- 
eracy and the Jahn- Teller nature of the relevant phonons, 
there is no contradiction between a phonon-mcdiated su- 
perconductivity and the relevance of electronic correla- 
tions, and the two interactions turn out to cooperate in 
providing relatively high critical temperature Q- 

In an extremely broad sense these materials (cuprates 
and fullerides), as well as many others that we did not 
talk about, raise the same conceptual problem, namely 
the investigation of systems in which both e-e interac- 
tion and e-ph coupling are non negligible and the physics 
can be explained only taking both into account. On the 
other hand the same phenomenology suggests that this 
competition may result in completely different physics 
according to specific aspects of the materials. In gen- 
eral we can expect different behaviors because of: (i) 
Different parameters within the same model (e.g., which 
is the largest scale between electron-phonon interaction 
and electron-electron repulsion); (ii) Different form for 
the interaction term, or more generally, different models. 

If we focus on point (i) we immediately realize that 
even the simplest models must contain several relevant 
physical parameters. As we will discuss in the following, 
we have to deal at least with the electronic bandwidth, 
the Coulomb repulsion, the strength of electron-phonon 
interaction, the phononic frequency and the chemical po- 
tential that controls the band filling. This determines a 



multi-dimensional phase diagram, which can hardly be 
understood in its entirety within a single analysis and it 
is expected to present several different regimes. There- 
fore, even if we choose one given simple model, it may be 
useful to focus on a given physical regime, which essen- 
tially implies to select a hierarchy between the different 
energy scales, or to fix (or neglect) some of them. 

Our choice, is to focus on the "strongly correlated" 
metallic phases, i.e., on system in which the Coulomb re- 
pulsion is the largest energy scale. The polar star of this 
work is the understanding of the fate of electron-phonon 
interaction in systems that are dominated by electron- 
electron interactions such as the cuprates. Nonetheless, 
our discussion will also follow some detours, which will 
help us to build a more comprehensive picture of the com- 
petition beteen the two interactions. One of these detours 
will touch point (ii) addressing the role of the phonon 
symmetry in its interplay with correlations. This point 
is crucial for the understanding of the sinergy between 
e-ph interaction and e-e correlation in the fullerenes. 

The manuscript is organized as follows: In Sec. II we 
introduce the Hubbard-Holstein model. Sec. Ill is de- 
voted to a Fermi-liquid analysis of the effects of correla- 
tions on electron-phonon interactions and to a mean-field 
solution of the Hubbard-Holstein model within the slave- 
boson formalism. Sec. IV presents a non-perturbative 
Dynamical Mean-Field Theory study of the Hubbard- 
Holstein model. Sec. V is dedicated to the charge 
instabilities of the model; Sec. VI briefly compares the 
Hubbard-Holstein model with a three-band model with 
Jahn- Teller interactions introduced for the fullerenes. 
Sec. VII presents our conclusions. 

THE MODEL 

The simplest model of a strongly correlated electron 
system coupled to the lattice is given by the single-band 
Hubbard-Holstein (HH) model 

- fi ^ n i° + U E n 't n 4 
icr i 

+ ujo^ a \ a i + 9^2 ( a l + a i) ( n i<y - ("iff)) 1 (!) 

i i,a 

where indicate nearest-neighbor sites. ni a = c| ff Ci ff 
is the local electron density, which is coupled via g to 
the field Oj of a dispersionless phonon. This model is 
clearly a very idealized description of a solid: both the 
Coulomb repulsion and the electron-phonon interaction 
are assumed to be purely local, only one electron band is 
considered, only the electron density couples to phonons. 
The relevant physical parameters are the strength of the 
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coulomb repulsion U, the dimensionless electron-phonon 
coupling A = 2g 2 N Q /oj (here N is the electronic density 
of states at the Fermi level), the adiabatic ratio cu /W 
(W being the electronic bandwidth), and the value of the 
chemical potential and the details of the bandstructure 
(inclusion of next-neighbor hopping). The dimensional- 
ity of the system also plays a major role. This multi- 
dimensional parameter space leads to an extremely rich 
physics, and the number of paper devoted to this simple 
model is countless. Various approaches have been con- 
sidered to solve the HH model in the presence of strong 
correlations. Besides numerical techniques like quantum 
Monte Carlo, exact diagonalization. |1 31-41 51] Dy- 

namical Mean Field Theory (DMFT), flil-1251]. Densitv- 
Matrix Renormalization Group [26f . (semi) analytical ap- 
proaches like slave bosons (SB) I27H31B . large-N expan- 
sions (32l . [33| . variational approaches |34j including a mod- 
ified Gutz wilier scheme [35l437j . have been useful to elu- 
cidate the renormalization of the e-ph coupling in the 
presence of (strong) correlations. 

According to our choice, in this work we will mainly 
consider on the large- U regime, and on regions of pa- 
rameters in which the system is a strongly correlated 
metal. Even if our focus will be the strongly corre- 
lated HH model, we will discuss its results in comparison 
with some related models, like the Hubbard-tJ model, 
the three-band Hubbard model for the cuprates and a 
three-orbital Hubbard model with Jahn- Teller interac- 
tions for the doped fullerides. Our investigation will be 
mainly dedicated to the effects of e-ph interaction on the 
self-energy and the quasiparticle renormalization factor z 
starting from a strongly correlated regime, in which the 
Coulomb interactions is strong enough to determine the 
physics of the problem, but the system is not a Mott in- 
sulator. Since the study of superconductivity in DMFT 
(our main tool of investigation) is limited to s-wave sym- 
metry we will not discuss superconductivity in the HH 
model which is expected to be d-wave if correlation dom- 
inates. On the other hand we will study s-wave super- 
conductivity in a three-orbital model for fullerides which 
emphasizes the role of the symmetry of e-ph interaction 
in presence of correlations. 



EFFECT OF ELECTRON-PHONON 
INTERACTION IN A CORRELATED METAL 

Fermi-liquid analysis 

In this section we begin our analysis of the properties 
of e-ph interaction in a correlated metal within a Landau 
Fermi-liquid (FL) picture [H, H^. Within this approach, 
the correlated metal is described as a collection of quasi- 
particles with an effective mass to* instead of the physical 
electron mass to. In the presence of strong e-e correla- 
tions the motion of the carriers is naturally obstacled by 



the interactions, which is reflected in a large ratio be- 
tween the effective mass and the bare mass to* /to 3> 1 
and in a loss of low-energy spectral weight described by 
a small quasiparticle renormalization factor z e . The for- 
mer reflects in an enhanced quasiparticle density of states 
N* = m* /mNo (No being the bare density of states) and 
the latter renormalizes the quasiparticle interactions. 

In order to characterize the fate of the e-ph interaction 
in a similar correlated metal we need to consider also 
the vertex corrections introduced by e-e interactions, for 
which no Migdal theorem can be invoked. 

We can gain a first insight on the way in which the e-ph 
interaction behaves in the presence of strong correlations 
by considering the effective dimensionless e-e interaction 
mediated by the exhange of a single-phonon 



N*Tf f M,w) = N*g'z e 'h e \q,u>)- 



eff 



(2) 



by assuming that all the correlation effects are local, or 
equivalently that the self-energy is independent on mo- 
mentum, the effective mass is related to z e by z e = 
m/m*. In Eq. @ A e includes the vertex corrections 
which renormalizes the e-ph vertex, and the last factor is 
the free phonon propagator. In order to focus on the ef- 
fect of correlations on the phonon-mediated interaction, 
in both z e and A e we include only processes due to e-e 
interactions, as reminded by the index "e" that we at- 
tached to them. 

Within a Landau Fermi-liquid picture we can use the 
Ward identities that connect the vertex corrections A 
with the wavefunction renormalization z. In the small 
frequency (w) and transferred momentum (q) regimes 
these identities have two distinct forms depending on the 
order of the u> — > and q — > limits. In the case of the 
charge-density vertex, which is relevant for our Holstein 
coupling, we have 



z e K e (uj -> 0, q = 0) = 1 
z e A e (ui = Q,q -> 0) = 



1 



(3) 



where Fq is the symmetric Landau parameter. ([3]) are 
exact Ward identities, which are satisfied irrespective of 
the details of the e-e interactions and show the drastic 
difference between the dynamic [(w — > 0, q = 0)] and 
static [(ui = 0, q — >• 0)] limits. 

Plugging these results into Eq. ^ we obtain, in the 
two limits considered above 



JVT# f (w-X), 9 = 0) 



N*T$ f (u = 0,q-K>) = 



2g 2 N* 
2g 2 N* 



2g 2 K e2 



(4) 



(5) 
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where n e = N* /[I + (-Fo)e] is the charge compressibil- 
ity in the absence of e-ph coupling. The difference be- 
tween the dynamic and static case can be dramatic in 
the case of a Fermi liquid with a large mass enhance- 
ment m* /to 3> 1, and small compressibility renormaliza- 
tion (k 6 ~ Nq), This requires {Fq) e to be much larger 
than one and proportional to the quasiparticle density of 
states N* = (m*/m)N > N . Eq. © then leads to 

iv*r^>^o, q = o) = -x(Tl) (6 ) 



N*T p J} f (id = 0,q^0) = -A 



(7) 



so that the effective one-phonon mediated e-e interaction 
is large (~ m*/m) in the dynamic limit and small (~ 
to/to*) in the static one. We emphasize that the key 
condition for the equalities (6,7) to hold is that x e "C 
77i* /to. Therefore they arc verified also for a modest mass 
enhancement as long as the compressibility remains much 
smaller than it. 

The strong uj — q dependence in Eqs.(6,7) has been 
demonstrated on really general grounds only in the small- 
q and small-w limits, whereas the case of finite g's and 
id's needs (approximate) analyses of specific models. The 
cases of the single-band and of the three-band Hubbard 
models with infinite local repulsion have been extensively 
considered in the literature as prototypical models of 
strong correlations. In this framework the issue of e-ph 
couplin g ha s been considered by means of the Holstein 
H HI, H3| or (less frequently) of the so-called 

Su-Schricffcr-Hccger couplingfwhere phonons couple to 
the electron hopping term) [40(. 

Results for these models show that the product zA re- 
mains of order one in the dynamical regime as long as 
the momentum and frequency lie outside the particle- 
hole continuum, while it is strongly suppressed (as in the 
static limit) inside the particle-hole continuum, where im- 
portant screening processes take place. Moreover, strong 
correlations provide further suppression of the e-ph cou- 
pling when, within the static limit id = 0, the transferred 
momentum is increased 32j . These additional screen- 



ing channels depend on the details of the electronic band 
structure determining particle-hole screening processes. 

The general Fermi liquid discussion and the specific 
analysis of models with strong correlations generically 
demonstrates the relevant role of dynamics in the screen- 
ing effects that e-e correlations induce on the e-ph cou- 
pling. This strong dependence of the e-ph vertex on 
momentum and frequency (and on their ratio) makes 
the effects of the e-ph coupling rather subtle, since dif- 
ferent physical quantities, involving different dynamical 
regimes, may display more or less suppressed e-ph ef- 
fects. In particular the e-ph coupling (and the e-e inter- 
action mediated by phonons) will be depressed by strong 
e-e interactions whenever small energy and large momen- 
tum transfer are involved (like, e.g., in transport). This 



suppression may be substantial, for instance, in the low- 
doping region of the superconducting cuprates, where 
e-e correlations are strong due to the relative proxim- 
ity to a correlation-induced insulating phase. On the 
other hand different physical processes involving dynam- 
ical processes could experience a more pronounced e-ph 
coupling. Specific calculations carried out in a single- 
band Hubbard-Holstcin model within a large-N treat- 
ment of the e-e correlations 32 , [33| demonstrated that 
the Eliashberg spectral function a 2 F(u>) determining su- 
perconductivity is much less reduced than the analogous 
quantity a 2 F tr (id) determining transport properties. As 
we will discuss in the following, this different renormal- 
ization will found a counterpart in nonperturbative dy- 
namical mean-field theory calculations. 

Even if our focus is on the Hubbard-Holstein model, it 
can be useful to recall that in the case of phonons coupled 
to the electron current, Ward identities similar to those 
of Eq. © can be derived [H| 



zA a (uj -+0,q = 0) 
zA a (td — 0, q -> 0) 



(8) 



th component of the electronic vec- 



and J qa are instead the a-th com- 



where A Q is the a 
tor vertex part. v q 
ponents of the quasiparticle velocity and of the current 
respectively. Similarly to the case of a coupling with the 
density, the correlations suppress much more strongly the 
e-ph coupling in the static limit, while they affect little 
the coupling in the dynamical case (remember that the 
current is weakly touched by interactions and it even re- 
mains constant in translationally invariant systems). 



The Hubbard-Holstein model: mean-field 
slave-boson approaches 

A natural tool to adreess the screening effects be- 
yond the small-q and small-w regime are mean-field ap- 
proaches based on a slave-boson representation of the 
Hilbert space. These methods, although approximate, 
capture the main physical ingredients of the problem with 
a description of a Fermi liquid of quasiparticles coupled 
by a residual interaction. This is why we briefly sum- 
marize here some results obtained [28| in the simplest 
formulation of the slave-boson large-N approach to the 
infinite-U Hubbard model. Here one can reach a semi- 
quantitative understanding of the effects of correlations 
on the e-ph coupling, which are in substantial agreement 
with the results of more sophisticated approaches. 

In this discussion we will consider the infinite-repulsion 
limit, which simplifies the formalism. In this limit we 
have a sharp contraint of no double occupancy on each 
lattice site n i<y — 1- The standard slave-boson tech- 
nique implements the constraint (41H45j | by perform- 



ing the usual substitution c ity 



->• 



b\fia, 



where the fermionic fi a operators represent quasiparti- 
cles, while the bosonic field &j labels the state of a site 
with no fermions in it. This formulation of slave bosons 
is here used in connection with a large- N expansion 42 1 
in order to introduce a small parameter allowing for a 
systematic perturbative expansion without any assump- 
tion on the smallness of any physical quantity. Within 
the large- TV scheme, the spin index runs from 1 to N and 



b]h = f . A 



the constraint assumes the form J2 a C L C 
suitable rescaling of the hopping amplitudes Uj — > t y /N 
must, in this model, be joined by the similar rescaling 
of the e-ph coupling g — > g/VN in order to compen- 
sate for the presence of N fermionic degrees of free- 
dom. It is beyond the scope of this review to report 
the technical details of this technique, which have been 
extensively presented in previous works [281 ] - Here we 
simply remind that the leading order in the large-N ex- 
pansion provides a mean-field description of the infinitc- 
U Hubbard model with uniform constant values of the 
bosonic field bi = bo and of the Lagrange multiplier field 
Xi = Xq implementing (on the average) the no-double- 
occupancy constraint. This gives rise to an insulating 
phase at half-filling (n = 1, doping 5 = 0) and a Fermi 
liquid metallic phase at finite doping with small quasi- 
particle residue z = b\ = b and large mass m* = 1/z. 
The treatment of the fluctuations at the next leading 
order beyond mean-field introduces residual interactions 
between the quasiparticles and allows to determine the 
scattering amplitudes in the particle-particle (Cooper) 
channel r c (k,k';q,uj) and in the particle-hole channel 
r(fc, k'\ q,uS). Taking the small-(q,oj) limits in this latter 
quantity also allows to determine the Landau parameter 



Fa = N*T U = UN*e kf 



N* 



(9) 



Here £k F is the electron dispersion calculated at the 
Fermi energy. As we already noted, Ffi enters in the 
FL expression of the compressibility. When Fq < — 1 the 
thermodynamic stability condition k > is violated and 
the system undergoes phase separation. We shall discuss 
this issue in Sec. VI. We anticipate that, in the phase sep- 
arated region, long-range Coulombic forces play a crucial 
role. It is indeed natural that a long-range interaction 
frustrates the formation of charge-rich regions. The out- 
come of this competition is, as we shall see, a shift of the 
charge instability to finite momenta promoting the for- 
mation of a charge-density-wave phase. In this section 
we will consider parameters far from the phase separa- 
tion instability. Nonetheless, in light of the central role 
played by long-range interactions in the phase separated 
case, we will also comment about their effect in these 
stable regions of parameters. 




FIG. 1: Leading-order in 1/N diagrammatic structure of the 
effective e-ph vertex dressed by electronic processes only: the 
dashed line is the slave-boson propagator only involving b and 
A bosons (pure e-e interaction), the solid dot is the dressed 
e-ph vertex, the open dots are the bare e-ph vertices and the 
grey dots are the quasiparticle-slave-boson vertices. 



Static properties 

In the SB large-N approach, the fluctuations of the 
bosonic fields mediate the residual interaction between 
the quasiparticles. If one only considers the fluctuations 
of the b and A fields, then only the effects of the electronic 
Hubbard repulsion are described and one can accordingly 
discuss the effect of the pure electronic screening preo- 
cesses on the e-ph vertex. For the specific model that we 
here briefly discuss, Fig. [ljreports the Feynman-diagram 
representation of the electronic processes (schematized 
by the dashed line of the bosonic propagators) dressing 
the bare (empty dot) e-ph vertex. The ratio between the 
resulting screened e-ph vertex and the bare e-ph coupling 
g is reported in Fig. [2] both in the absence (dotted curve) 
and in the presence (solid curve) of a long-range Coulomb 
repulsion. In this latter case the strength is chosen to 
produce a repulsion of about O.IeV between electrons in 
nearest-neigbor cells. In both cases there is a strong re- 
duction of the static e-ph coupling. In the short-range 
case the residual coupling is strongest at small momenta 
and it decreases as the momentum. 

On the other hand, as it is natural, the long-range 
potential screens out the long-range charge fluctuations 
thereby driving to zero the e-ph coupling at low mo- 
menta. 

These findings can be obtained and confirmed with 
several different approaches. In particular, they repro- 
duce exactly the results of large-N calculations with 



Hubbard projectors instead of SB's [32], |33|, and are 
in good agreement with calculations based on the flow- 
equation method[46| and recent Gutzwiller+RPA calcu- 
lations at finite (but large) u[37j]. They also agree quali- 
tatively QMC analysis [11( (although some differences are 
present, which can be attributed to the fact that QMC 
calculations are performed at finite temperature and at 
finite Matsubara frequencies) (cf. Figs. 11 and 12 of Ref. 
33). 

All these results show that e-ph scattering at large mo- 
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FIG. 2: Static effective e-ph vertex in units of the bare e- 
ph coupling g as a function of the transferred momentum (in 
units of the inverse lattice spacing l/o) in the (1,0) direction. 
The vertex is calculated at leading order in 1/JV for a HH 
model with t = 0.5eV, t' = -(l/6)t, ujo = 0.04eV and doping 
x = 0.205. The dotted line is in the absence of long-range 
Coulomb forces (Vc = 0); the solid line is in the presence of 
long-range Coulomb forces with Vc = 0.55eV (adapted from 
Ref. [H) 

menta is typically weaker than scattering at low momenta 
in the presence of short range forces only (which is the 
case of metallic phases far from phase separation charge 
instabilities, where, on the contrary, long-range interc- 
tions start to play a crucial role). This can be of obvious 
relevance when the relative importance of e-ph couplings 
between the quasiparticles and specific phononic modes 
is considered [471 ]. Indeed it might well happen that in 
the presence of strong correlations modes that would be 
strongly coupled, but would exchange preferably large 
momenta, are more severly screened than other modes 
exchanging lower momenta. All this surely deserves a 
specific analysis. 



Dynamical properties 
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FIG. 3: Effective e-ph vertex in units of the bare e-ph coupling 
as a function of the Matsubara frequency in the HH model at 
leading order in 1/N. The parameters are the same as in Fig. 
[2] Panel (a) is for a small transferred momentum q = (0.2, 0) 
(in units of inverse lattice spacing 1/a); (b) is for a sizable 
transferred momentum q = (2.0,0) (after Ref. [28l |1 

Clearly in the former case the e-ph coupling rapidly re- 
turns to its bare value, at Matsubara frequencies larger 
than ~ up|q| . The vp\q\ scale is instead much larger 
in panel (b), where the effective e-ph vertex stays small 
over a much broader frequency range. This fully parallels 
the low-frequency and low-momentum limits discussed in 
Sect. III.A. 

Again this result is not specific of this single-band 
model (or of its treatment) and it has been confirmed 
by the analysis of a three-band Hubbard model for the 
cuprate Cu02 planes with infinite repulsion on copper 
orbitals. In this case, one can even notice in the small- 
momentum case an enhancement the effective e-ph cou- 
pling above its bare value. This is an overscreening effect 
due to interband processes. In any case, again one finds 
that low-momentum processes generically lead to larger 
e-ph couplings. 



The previous subsection was focusd on the screening of 
the static e-ph coupling by electronic processes. However 
the Fermi-liquid analysis carried out in Sec. „ pointed out 
the great difference between screening processes in the 
static and in the dinamical regimes. While the Fermi 
liquid analysis was only able to provide definite state- 
ment in the small-^, u>) regime, within the SB large-N 
approach we can investigate the role of dynamics in the 
screening processes beyond this limit. 

Fig. [3] displays the behavior of the effective e-ph ver- 
tex (again normalized to the bare g) as a function of 
Matsubara frequencies for two distinct momenta in the 
(1,0) direction for our HH model in the infinite-?/ limit. 
In panel (a) a small momentum q = (0.2, 0) (again unit 
lattice spacing is used here) is reported, while panel (b) 
shows the behavior at a larger momentum q = (2.0,0). 



STRONG-COUPLING REGIME AND POLARON 
FORMATION IN A CORRELATED METAL 

Proximity to a paramagnetic Mott insulator 

In this section we extend our analysis behind the mean- 
field level and discuss the fate of the e-ph interaction in 
a correlated metal under the sole assumption that the 
physics is governed by the Hubbard repulsion, without 
assuming a weak e-ph coupling, and/or any approxima- 
tion as far as the adiabatic ratio is concerned. This 
means that we need a theoretical approach able to treat 
severeal energy scales symultaneously, without assuming 
that any of them is negligible or perturbative. A natural 
candidate for this purpose is the Dynamical Mean-Field 
Theory (DMFT)[48|. which treats all local interaction 
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terms (such as both the Hubbard and the Holstein cou- 
plings) and the hopping term on the same footing and it 
is equally suitable to treat any parameter regime. The 
central approximation beyond DMFT is the locality of 
the self-energy (both the electronic and the phonic con- 
tributions), a condition which becomes exact when the 
coordination number becomes large. The original lat- 
tice enters in the calculation only through the density of 
states, which we always choose to be a semicircular one 
of width W. 

DMFT allowed to obtain a complete characterization 
of the Mott-Hubbard transition in the pure Hubbard 
model, and the emerging physical picture is able to ex- 
plain several properties of correlated oxides. While we 
refer to original papers 48| for details, we recall here some 
aspects which are relevant to our discussion. 

We first consider the situation in which the number of 
electrons equals the number of lattice sites (the so-called 
"half-filling" situation) and, just to focus on pure corre- 
lation effects, we consider a paramagnetic phase. In this 
regime, for large repulsion, the ground state of the HH 
model can become a "Mott" insulator, in which the elec- 
trons are localized because the electron motion is ener- 
getically unfavorable. Starting from the uncorrelated sys- 
tems and increasing the correlation strength, the spectral 
weight is transfered from low to high frequency. In this 
process the spectral function evolves from a single band 
to a three-feature structure in which a renormalized band 
survives around the chemical potential, while precursors 
of the Hubbard bands develop around oj ~ ±[7/2. As the 
correlation is further increased up to a value U = U C 2 the 
quasiparticle band disappears, and the system becomes a 
Mott insulator with a preformed Hubbard gap. The key 
parameter that controls this Mott-Hubbard transition 
from a metal to an insulator is the quasiparticle weight z, 
which is directly computed from the self-energy (which is 
a natural outcome of a DMFT calculation), z measures 
the width and the total spectral weight of the low-energy 
quasiparticle peak, and its vanishing pinpoints the Mott 
transition. As soon as the system is doped away from 
half-filling, the quasiparticle peak moves towards one of 
the bands, but it is still well defined for a sizeable dop- 
ing region (dependent on the value of U and on details 
of the bandstructure) [IH . Due to the momentum inde- 
pendence of the self-energy, also in DMFT we have that 
z = m/m*, so that the Mott transition is associated with 
a divergent effective mass of the renormalized carriers. 

In the following we discuss the effect of a non- 
perturbative electron-phonon coupling in the strongly 
correlated metallic solution, i.e., a system in which a 
quasiparticle peak at the Fermi level is separated from 
the Hubbard bands realizing a separation of energy scale. 
Starting from this situation, in which the quasiparticle 
bandwidth is z times the bare width W, the effect of the 
e-ph coupling is far from trivial. Indeed there are two 
main effects associated to e-ph interaction: 



The e-ph interaction can introduce a further quasi- 
particle renormalization, associated to the increase 
of the quasiparticle effective mass, which may even- 
tually lead to polaronic effects for very strong cou- 
pling. This effect leads to a decrease of z. In weak 
coupling we have l/z~l + ^Xat half-filling and a 
semicircular density of states of half width ly.jlij] 

The e-ph interaction mediates an attractive 
density-density interaction (the density-density 
form is specific for a Holstein coupling), which di- 
rectly contrasts the Hubbard repulsion. If we in- 
tegrate out the phonon degrees of freedom, the 
fermions interact through a dynamical (retarded) 
interaction. 

2g 2 u) 



U eff {u) = U- 



(10) 



In the antiadiabatic regime the frequency depen- 
dence of the second term can be neglected and over- 
all interaction is a static term J7 stQt = U — XW, in 
which the e-ph interaction reduces the strength of 
the Hubbard term. When the phonon frequency 
becomes finite the interaction is retarded, but we 
still expect a similar effect. Such a decrease of the 
effective repulsion is expected to make the system 
less correlated and increase z (assuming that the 
repulsion still overcomes the attraction) 

The balance between this two effects is not generic 
and it depends on the adiabatic ratio and on the pre- 
cise value of the interactions. Yet, important conclusions 
can be drawn in the correlated regime, in which, also in 
the presence of e-ph interaction, the separation of energy 
scales determined by correlations survives. In this re- 
gard, it is important to recall that, within DMFT, the 
quasiparticle weight is associated to a Kondo resonance 
of an Anderson-Holstein impurity model. Assuming that 
the Hubbard U is the largest scale of the problem, the 
Kondo coupling can be calculated in terms of virtual pro- 
cesses acting in the subspace in which the impurity is 
singly occupied obtaining an effective Hamiltonian for 
spin fluctuations [H(| ■ The result is given by 



Jk(X) = Jk(0) 



oo |^o|e ff / tJo ^ a_at ^| 
^ 'l-2g 2 /uj U 

m=0 a ' 



2muj /U 



(11) 



where the Kondo coupling in the absence of phonon is 



given by J K (0) = 16V 2 /U, 



is the state with m 



phonons, and V is the hybridization between the impu- 
rity and the bath. After some algebra, and introducing 



U eff = U-r)\W, 
we can write, for small XW/U 

XW\ 16V 2 



(12) 



J K (X)~J K {0)[l+ri 



U 



16V 2 



U eff U -rjXW 



(13) 
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FIG. 4: Effect of the electron-phonon interaction on the quasi- 
particle weight z in the presence of electron-electron interac- 
tion of different strengths. The quantity plotted as a function 
of A is the ratio between the full z and the value of the Hub- 
bard model without e-ph coupling (z(0)) (after Sangiovanni 
et al. HH). 



with 



V = 



2u) /U 
1 + 2oj /U 



(14) 



This result would imply that the complicated interplay 
between the static Coulomb repulsion and the retarded 
e-ph coupling may be effectively described by an effec- 
tive purely electronic Hubbard model with a reduced re- 
pulsion. We notice in passing that this calculation for 
the Anderson impurity model is analogous to the evalu- 
ation of the effect of phonons on the superexchange cou- 
pling of Ref. 51]. Interestingly the phonon dynamics 
only enters thtough the ratio ujo/U, which can be consid- 
ered as typically small because U is by choice the largest 
scale of the problem, and uiq is smaller than the hopping 
scale. In the relevant regime of small coo/U, we have 
U eff = U-4g 2 /U 2 . 

These results can be tested through a DMFT solu- 
tion of the Hubbard-Holstein model at half-filling, where 
the separation of energy scales is more solid. In Fig. |4] 
we show the ratio z(X)/z(0), in order to emphasize the 
phonon contribution to the quasiparticle weight. It is ap- 
parent that the value of U determines different regimes. 
For small U the e-ph interaction reduces z (i.e., increases 
the effective mass) as expected in weakly interacting sys- 
tems. Increasing U we approach an opposite behavior 
in which the e-ph interaction makes the quasiparticles 
lighter [even if they are obviously heavier than free parti- 
cles because of the stronger renormalization determined 
by correlation, which is hidden in z(0)]. 

This behavior confirms that, when the e-e correlation 
dominates, the leading effect of the e-ph interaction is 
a reduction of the effective U, resulting in an increased 
quasiparticle mobility. We can now test the above predic- 
tion of an effective static repulsion including the effects 
of e-ph interaction as far as the low-energy physics is 



FIG. 5: Effective static electron-electron interaction for the 
low-energy properties of the Hubbard-Holstein model. The 
picture shows the coefficient r\ in U e f / = U — r/XW as a func- 
tion of wq/U (after Sangiovanni et al. 
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FIG. 6: Comparison between the DMFT spectra for the 
Hubbard-Holstein model and the effective purely electronic 
Hubbard model with U = U e ff- The low-energy part of the 
spectra is perfectly reproduced by the effective model, while 
the high-energy features present phononic signatures that can 
not be plugged in an effective static interaction (after Sangio- 
vanni et al. \2M). 



concerned. 

Assuming a form J7 e // = U — r]XW for such an effec- 
tive interaction, we obtained 77 for several values of U 
and wo simply determining the value of U which gives, 
for a pure Hubbard model, the same z we obtain for the 
Hubbard-Holstein model. The results are summarized in 
Fig. [5l and confirm that rj is essentially a function of 
loq/U which starts off linear at small values of the ar- 
gument before bending for larger values. The numerical 
value beautifully fits the value given above on the basis 
of the analogy with Kondo effect. 

On the basis of our knowledge about the DMFT re- 
sults for the pure Hubbard model, the quasiparticle 
weight completely characterizes the low-energy quasipar- 
ticle peak. Therefore our analysis should imply that the 
low-energy part of the spectral function of our Hubbard- 
Holstein model can be described by means of the effective 
Hubbard model that we introduced. This is strikingly 
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confirmed by a direct comparison, as shown in Fig. [6] 
Here we show some examples of spectra for the Hubbard- 
Holstein model compared with the associated effective 
Hubbard model. Besides the spectacular confirmation of 
the validity of the effective model for the low-energy part 
of the spectrum, another important information emerges 
from the picture: The high-energy part of the spectrum is 
instead affected by phonons in a more "dynamical" way, 
meaning that the high-energy Hubbard bands acquire a 
modulation in frequency which can be related to phonon 
satellites, which are completely absent in the low-energy 
part, where nothing happens at the characteristic phonon 
frequency. 

The picture that emerges from the previous arguments 
can be summarized as follows: Quasiparticle motion 
arises from virtual processes in which doubly occupied 
sites are created. Obviously, these processes are not so 
frequent, since the energy scale involved is large, but they 
are extremely rapid (the associated time scale is cx 1/U), 
and consequently are poorly affected by phonon excita- 
tions with a characteristic time scale 1/ujq » 1/U. When 
the phonon frequency is small with respect to U, the 
phonon degrees are frozen during the virtual excitation 
processes. Therefore, despite the overall electron motion 
is quite slow due to the small number of virtual processes 
(which is reflected by the large effective mass), the e-ph 
interaction has no major effect except for a slight reduc- 
tion of the total static repulsion. 

According to what we have just described, we can con- 
clude that strong correlations reduce the effect of the 
e-ph interaction on the low-energy properties, associated 
to quasiparticle propagation, while the high-energy prop- 
erties present more standard phonon signatures, such as 
the satellites at the phonon-frequency scale. 

It has to be underlined that DMFT is not able to intro- 
duce momentum-dependent corrections to the electronic 
properties. The above analysis therefore shows indeed 
that the "standard" electron-phonon interaction is heav- 
ily screened (and it actually loses its dynamical nature) 
when the low-energy quasiparticle properties are consid- 
ered. Only "non-standard" effects, such as the prevalence 
of forward scattering that we discussed in Sec. „can sur- 
vive at low energy. 

The above scenario has been carried out at half-filling, 
where the presence of e-e correlations has its most strik- 
ing effects, both in terms of the phase diagram, as it can 
give rise to a Mott transition, and in terms of the sep- 
aration of energy scales, which is clearly sharper than 
for doped systems. Therefore, as soon as we dope the 
Hubbard-Holstein model, even for U > A the dominance 
of repulsive correlations is weaker and the interplay with 
e-ph coupling will be more subtle [22| . As a result, it 
is not possible to describe the effect of phonons on the 
highly correlated metal in terms of an effective static po- 
tential unless the system is very close to the antiadiabatic 
limit. Therefore, for small and intermediate phonon fre- 
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FIG. 7: Effective mass m* as a function of A for three different 
densities different from half-filling (n=0,5, 0.7, 0.9) in the 
uncorrelated (left) and strongly correlated (U = 2.5W , right) 
system (after Sangiovanni et al. [221 ]') 



quencies, we do not find situations in which increasing 
the electron-phonon coupling reduces the effective mass. 

Nonetheless, if we consider reasonably large values of 
U, the dominance of correlations will determine a reduce 
effectiveness of the e-ph coupling, and, e.g., polaron for- 
mation will be pushed to significantly larger values of A 
then for uncorrelated systems, as shown by the DMFT 
results of Fig[7] Here we plot to* /to as a function of A 
and we compare the uncorrelated system with the sys- 
tem with U/W— 2.5 (Here, since the density is different 
from half-filling, the system is always metallic even if 
U is larger than the critical value for the Mott transi- 
tion). While in the uncorrelated case m* grows exponen- 
tially when A approaches a critical value of order I for 
all densities, signaling polaron formation, the correlated 
system displays a significantly weaker growth of to* up 
to A ~ 0.5 ~ 0.75. 

The role of the adiabatic ratio is illustrated by Fig. 
[SJ where we report the renormalization of the linear co- 
efficient of the mass enhancement defined by the rela- 
tion m*(U)/m*(U, A) = 1 - rX. Here m*(U, A) is the 
effective mass in the presence of both electron-electron 
and electron-phonon interactions and m*(U) is the same 
quantity in the absence of coupling to the phonons. Here 
a negative r implies a standard increase of the effective 
mass due to phonons. The results (again for U/W=2.5) 
show that in all cases the coefficient is smaller than one, 
confirming that correlations reduce the effective e-ph cou- 
pling, and a strong (and nonmonotonic) dependence on 
the antiadiabatic ratio. Only for very large values of 
ujo/W r becomes positive reflecting that the "screening" 
physics we described above is effective. The dependence 
on the density naturally reflects that fillings closer to 
n = 1 display weaker phononic effects in the adiabatic 
regime of small frequencies and a more rapid evolution 
to an antiadiabatic regime in which r > 0. 

The scenario which emerges from DMFT calculations 
at finite U /W can be confirmed by a semi-analytical ap- 
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FIG. 8: Renormalization of the quasiparticle electron-phonon 
coupling as a function of ujq/W for three different densities 
and U/W = 2.5 (after Sangiovanni et al. [22T |) 




FIG. 9: Renormalization of the effective mass due to electron- 
phonon coupling for infinite U and ujo/W = 0.1 as a function 
of A for different doping levels. The inset shows the same data 
in a smaller range. Here A is twice as in the rest of the paper, 
(after Barone et al. [3^|) 



proach based on an extension of the Gutzwiller approach 
which treats phonon degrees of freedom on the same foot- 
ing of the electrons. For more details see Refs. [3|1 l3rjj ]. 
The results of Fig. |9]in the limit of infinite repulsion con- 
firm the scenario arising from Fig. in the correlated 
system the electron-phonon interaction is considerably 
reduced with respect to the free system, but the qualita- 
tive effect of A is an increase of the effective mass. 



Antiferromagnetic correlations 

The above analysis suggests that, as expected, strong 
e-e correlation essentially opposes to e-ph coupling, even- 
though "anomalous" signatures of e-ph coupling still sur- 
vive even in regions of parameters for which correlations 
prevail. Yet, these results are limited to the metallic 
paramagnetic state, in which no broken symmetry is al- 
lowed. At half-filling and for some finite doping region, 
strong correlations lead to an antiferromagnetically or- 
dered state, and it is expected that finite-range antifer- 



romagnetic correlations survive in a wider doping region. 
The relation between antiferromagnetism and e-ph inter- 
action is hinted by the experimental framework. Indi- 
cations of phonon signatures in high-T c superconductors 
are indeed particularly strong in the extremely under- 
doped region, where some kind of antiferromagnetic cor- 
relation is certainly present. For example, clear polaronic 
features are observed in the optical spectroscopy^ and 
ARPES[3] of underdoped materials. 

From a theoretical point of view, several investiga- 
tions indeed suggest that the e-ph interaction is par- 
ticularly effective for a hole in an antiferromagnetic 
background 53l 54 [, and for slightly doped t-J models. 
These results can be reconciled with the above findings 
for the nonmagnetic phase by simple arguments. 

Within the paramagnetic phase, the effect of increas- 
ing correlations is a strong reduction of the quasiparticle 
weight z associated with a divergent self-energy, which in 
turn strongly renormalizes the e-ph vertex, leading to the 
strong reduction of low-energy phononic signatures that 
we described above. Once antiferromagnetic correlations 
are allowed, the system can turn insulating even with a 
finite z and a non-divergent self-energy, hence the e-ph 
vertex is not severely screened. From a more physical 
point of view, the antiferromagnetic insulator allows for 
more charge fluctuations with respect to the pure Mott 
state at the same value of U. Therefore a Holstein cou- 
pling, which exploits precisely charge fluctuations to gain 
energy, is expected to be favoured by antiferromagnetic 
correlations. 

Direct DMFT calculations in the antiferromagnetic 
phase at half-filling confirm these expectations. In Fig. 
[TO]we show the quantity [z/z(0) — 1]/A for small A {z(0) 
being z in the absence of e-ph interaction) . This quan- 
tity measures the renormalizaton of the linear e-ph cou- 
pling induced by e-e correlations. The comparison be- 
tween the paramagnetic solution and the antiferromag- 
netic state confirms the above expectations. While this 
coefficient rapidly drops as a function of U in the param- 
agnetic state, the inclusion of antiferromagnetism leads 
to a much more robust e-ph coupling. We underline that, 
however, the e-ph interaction is still substantially reduced 
with respect to the non-interacting systems. 

The comparison with the uncorrelated system is shown 
in Fig. [TTJ where the evolution of z as a function of A is 
followed beyond the perturbative regime. In all cases z 
decreases monotonically until it vanishes signaling po- 
laron formation. Yet, the decrease is more rapid for 
the uncorrelated system, and the critical coupling is sig- 
nificantly enhanced in the correlated antiferromagnetic 
state. 
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FIG. 10: Effect of the Coulomb repulsion on the linear con- 
tribution in A to the quasiparticle renormalization. DMFT 
results for the antiferromagnetic phase (AF-DMFT) are com- 
pared with DMFT in the paramagnetic state (P-DMFT) and 
to a self-consistent Born Approximation (SCBA) for the t- 
J model for parameters corresponding to the red dots AF- 
DMFT. Antiferromagnetism allows for a sizeable e-ph cou- 
pling also in the present of very strong Coulomb repulsion, in 
constrast with paramagnetic results (after Sangiovanni et al. 
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FIG. 11: Quasiparticle renormalization as a function of A 
for different values of U and cjo ~ 0.01251F. The critical 
value for polaron formation is only moderately increased by 
correlations with respect to the non-interacting case (after 
Sangiovanni et al. [551]') 



PHONON MEDIATED CHARGE INSTABILITIES 

In the previous sections we focused on the metallic 
phases, without discussing the possible instabilities, ei- 
ther directly driven by the interaction terms, or favoured 
by the weakness of the correlated metallic state. One can 
indeed expect, on very general grounds, that the reduced 
kinetic energy characteristic of the strongly correlated 
metal can be easily overcome by different localizing ef- 
fects thereby destabilizing the metal in favor of ordered 
phases. 

As we did in Sec. „in the discussion of the e-e screening 
of the e-ph vertex, we can start our analysis in the FL 
framework before discussing in some more details model- 



specific results. 

In particular, we focus on charge instabilities, which 
have been invoked to play an important role in high-T c 
superconductors. The charge compressibility k — dn/dj-i, 
which in a FL theory reads 



N* 



1 + F? 



(15) 



is the key quantity that controls the stability of the 
charge degrees of freedom. A positive k is the stabil- 
ity condition. Therefore if Fq < —1 this condition is 
violated and the system undergoes a charge instability 
towards a phase-separated state. 

It is important to stress that the s of Sec. „about the re- 
duction of the static e-ph vertex do not necessarily imply 
that e-ph effects do not contribute to the charge com- 
pressibility in the presence of strong correlations. Indeed 
Fq = 2N*T UJ , where T w is the full dynamic effective scat- 
tering amplitude between the quasiparticles [38$ including 
both e-e and e-ph interactions. As opposed to its static 
counterpart, the dynamic amplitude which enters Fq is 
not depressed by the e-e vertex corrections. For example, 
at lowest order in g 2 and by performing an infinite-order 
RPA resummation of the e-ph screening processes, Fq 
reads 



F° = 2N* (TZ-g 2 



(16) 



where is the dynamic vertex function determined by 
e-e correlation processes only. This equation indicates 
that a sufficiently large bare g 2 can overcome the effective 
residual repulsion between the quasiparticles leading 
to a phase separation instability marked by the Pomer- 
anchuk condition F S < — 1 . It is worth pointing out that 
m*/m ^> 1 requires a large bare repulsion between the 
physical electrons (a large Hubbard U in our model) but 
this by no means requires a large residual interaction be- 
tween the heavy quasiparticles. Therefore even a small 
e-ph interaction can give rise to a phase separation insta- 
bility. 

Moreover, near the instability condition Fq = —1, the 
phonon contribution to the vertex becomes substantial 
and the e-ph interaction becomes relevant even in the 
static limit, at least at small g's. At large q the analysis 
of specific models shows that the e-e interaction mediated 
by phonons is instead suppressed also near the instability 
region. 

This is what indeed happens in the HH model treated 
within the SB-large-N method described In Sec. „ Within 
this approach it was first demonstrated that a metal with 
moderate e-ph coupling and strong e-e correlations could 



undergo a charge instability [28|, [56J . Specifically in the 
absence of long-range Coulombic forces the Pomeranchuk 
stability condition is violated. The doped HH model does 
not form a uniform phase and the system undergoes a 
phase separation between hole-rich regions and insulating 
half-filled regions. 
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FIG. 12: Phonon dispersion curves (a) in the instability 
(0.28,0.86) direction and (b) in the (1,1) direction, for t = 
0.5eV, t' = -(1/6)*, g = 0.240/V2eV, w = 0.04eV, and 
Vc = 0.55eV. The solid curves correspond to the critical dop- 
ing 5 C = 0.195, the dot-dashed and the dotted curves corre- 
spond to S = 0.21 and 5 — 0.3 respectively. 



In the presence of long-range Coulomb interactions the 
electrostatic cost of the charge-rich regions would become 
infinite, and the thermodynamic phase separation can 
not take place. The tendency towards charge segrega- 
tion would however survive and it finds its realization 
with a charge-density wave instability which emerges as 
a compromise between the charge segregation tendency 
and the homogeneizing effect of long-range interactions. 
This mechanism for charge ordering is the so-called "frus- 



trated phase separation" [57H60I ]. For the specific HH 



model described here it was found that for realistic val- 
ues of the e-ph coupling and of the long-range repulsion, 
frustrated phase separation gives rise to a second-order 
quantum phase transition (quantum critical point, QCP) 
around optimal doping (doping x = 0.19) [561 ] with the 
ordering and the periodicity influenced but not directly 
related to the structure of the Fermi surface. This in- 
stability arises instead from the energetic balance be- 
tween the tendency to phase separation and the frustrat- 
ing electrostatic cost of the long-range Coulomb interac- 
tion. Near this instability the phonon spectrum becomes 
highly anomalous. First of all the phonon acquires a 
strong coupling to the electronic degrees of freedom near 
the instability wavevector q c [which usually tend to oc- 
cur at q c « (±7r/2,0), (0, ±tt/2) for the relevant dopings 
and Fermi-surface shapes] . Near the instability wavevec- 
tor the phonon linewidth becomes therefore very broad 
and it even acquires a background of the order of the 
particle-hole continuum. At the same time the phonon 
dispersion softens and at the instability the frequency of 
this mixed phonon-electron mode vanishes. Remarkably, 
since the critical wavevector is not so large (typically of 
the order of 7r/2) the region where the phonon disper- 
sion becomes strongly anomalous is rather isotropic and 
substantial anomalies are present also in the (1,1) direc- 
tion upon approaching the critical doping of the QCP. 
Fig. [T^] reports the anomalous phonon dispersion found 
in the HH model in Ref. [28 1. 



Of course it is quite tempting to relate these anomalies, 
to the anomalies observed by inelastic neutron scattering 
6ll463{. An alternative possibility can also be proposed 
for the anomalied detected in underdoped cuprates: 
charge ordering can give rise to rather anisotropic nearly 
one-dimensional dynamical charge textures. In this case 
Kohn anomalies can be expected along the stripes at 
wavevectors of the order 2kp in the stripe direction 64 1 



JAHN-TELLER COUPLING IN THE 
FULLERENES 

All the above analysis has been carried out for the 
HH model. It has to be emphasized that some of the 
effects we discussed may be less general than what the 
simple form of the Hamiltonian may suggest. As we dis- 
cussed in details, the HH model is indeed characterized 
by two interaction terms which are both related to the 
charge degrees of freedom, and they indeed directly com- 
pete, as clearly shown by Eq. ([TU|). This direct competi- 
tion makes the two effects more exclusive than in general 
situations in which the e-ph coupling does not directly 
compete with the Coulomb repulsion. We can have two 
different ways to avoid the direct competition: a differ- 
ent functional form for the e-ph coupling in a single-band 
model, or a Jahn- Teller couping in a multiorbital model. 

The first situation can obviously have relevance for the 
cuprates, in which different phonon modes with specific 
smmetries may play a role, or for system dominated by 
the so-called Su-Schrieffer-Heeger coupling in which the 
phonons modulate the nearest-neighbor hopping. The 
second situation occurs instead in the fullerenes, where 
the relevant conduction band is a three-fold degenerate 
manifold of t\ u (p-like) orbitals, which couple with Jahn- 
Teller active local distortions of the fullerene molecule. It 
is precisely this kind of coupling which is expected to be 
responsible for supercondutivity in these compounds [tjBj]. 

A three-band model which includes a strong Coulomb 
repulsion, a Hund's rule splitting and a moderate Jahn- 
Teller e-ph coupling has been studied in several papers 0, 
66 - 6^| , reaching an a priori suprising concousion: The 
Jahn- Teller coupling is not harmed by large Coulomb re- 
pulsion, and the phonon-driven superconductivity can ac- 
tually be strongly enhanced in the proximity of the Mott 
transition, i.e., in the region in which the correlations are 
most effective. 

Here we do not discuss the physics of this model in 
details, since we are mainly interested in contrasting its 
behavior with the HH model. The key point is that the 
Jahn- Teller interaction does not touch the total charge 
on each molecule, as it couples with a combination of 
local spin and orbital degrees of freedom. As a conse- 
quence, even when correlations are sufficiently strong to 
suppress the electronic motion, the localized electrons 
can still interact within a single fullerene molecule via 
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the e-ph interaction. For example, if we consider the 
experimentally relevant situation of three electrons per 
fullerene molecule, as we approach the Mott state three 
electrons will remain stuck on each fullerene. Yet, the 
can still be in a high-spin state or low-spin state, and the 
energetic gain associated to the multiplet splitting will 
the the same as for a non-interacting molecule. There- 
fore, the e-ph driven interaction will be not renormalized 
by correlations, as opposed to the Holstein model. From 
a Fermi-liquid point of view, the lack of renormalization 
is determined by very large vertex corrections (divegent 
like 1/z as the Mott transition is approached) that com- 
pensate the z factors. 

As a matter of fact, the effective interaction be- 
tween quasiparticles obtained in a nonperturbative 
DMFT study of the model corresponds to a severely 
screened Hubbard repulsion plus an essentially un- 
screened phonon-driven attraction that can be param- 
eterized as 



(17) 



This simple equation shows that, even if U is chosen to 
be significantly larger than J, when the Mott transition 
is approached (i.e., z — > 0)[(36j], the attraction will eventu- 
ally prevail. Moreover, in this regime the quasiparticles 
are quite heavy, and their large effective density of states 
can lead to an enhancement of the effective dimension- 
less coupling, which is expected to reflect in an enhanced 
critical temperature. 

This is precisely what is found by explicitly solving 
the superconducting phase of our three-band model in 
DMFT. If we follow the evolution of the superconduct- 
ing order parameter as a function of U for a fixed small 
J we first have a standard BCS-like region when U is 
so small that the bare attraction simply overcomes the 
bare repulsion. Then superconductivity disappears be- 
cause U is large enough to kill the attraction, but z is 
still close to 1. Further increasing U we approach the 
Mott transition, z decreases and it strongly renormal- 
izes down the effective repulsion. Eventually the effec- 
tive interaction becomes attractive and superconductiv- 
ity re-emerges, with an order parameter that follows a 
bell-shaped curve before vanishing at the Mott transition 
point. This strongly correlated superconducting pocket 
displays a maximum critical temperature which exceeds 
the weak-coupling BCS value. In other words, phonon- 
driven superconductivity is actually enhanced by strong 
correlations 0- 

A full DMFT solution of the model has allowed both to 
predict the experimental observation only later provided 
in Refs. @, like the dome-behavior of the critical tem- 
perature as a function of doping, and the first-order tran- 
sition to a spin-1/2 antiferromagnet when pressure is re- 
duced to recover the ambient phase of A15 CS3C60 , and to 
further characterize the properties of strongly correlated 



superconductors. For example, we predict a pseudogap 
in the photoemission spectra, and a kinetic-energy driven 
superconductivity for the most expanded compounds frjjj]. 

In the context of this review, our solution is a clear 
example of the crucial role of the phonon symmetry In 
our multiband model it is possible to consider phonons 
which are by symmetry unharmed by correlations, as op- 
posed to the Holstein model. The result is confirmed 
by investigations of simplified two-orbital models which 



share the same properties 69, 70] When we go back to the 



cuprates, and to single-band models, our findings sug- 
gest that phonon modes which are coupled to operators 
which are not proportional to the charge can indeed sur- 
vive much better in a strongly correlated environment, 
in analogy with the findings of mean-field methods. 



CONCLUSIONS 

The focus of the present topical review is on the effects 
of strong electron-electron correlations on the electron- 
phonon coupling. We mostly considered the Hubbard- 
Holstein model, where both e-e and e-ph interactions lo- 
cally couple to electron density fluctuations. In this case 
the competition between these interactions is quite effec- 
tive, particularly in the proximity of the Mott-Hubbard 
transition. Indeed, since the e-e Hubbard repulsion 
makes the density fluctuations stiffer, the Holstein e-ph 
coupling g is generically suppressed. Along this paper 
we made more specific this point by reviewing the spe- 
cific momentum and frequency dependencies of this sup- 
pression. A first observation is that the suppression is 
strong whenever the quasiparticle residuum z is small. 
On the contrary we observed that in the antiferromag- 
netic phase, where the suppression of double occupancy 
imposed by the large Hubbard repulsion U is due to the 
spin ordering and does not entails a small z, the e-ph 
coupling is only weakly suppressed. This easily explains 
why polaronic features are present and clearly visible in 
weakly doped antiferromagnetic cuprates. The weakness 
of this suppression is likely to persist even in the metal- 
lic paramagnetic regime, if substantial residual antiferro- 
magnetic correlations persist on a local basis. The oppo- 
site case of a strong suppression of the e-ph coupling Re- 
curring for small z (or more precisely when n e /N* <C 1), 
needs further specification. In particular we find that g 
is more or less suppressed depending on the dynamical 
regime: for a small ratio between the transferred fre- 
quency uj and the transferred momentum vfQ the e-ph 
coupling is strongly reduced, while in the opposite limit 
ui/vFq_^> 1 no suppression is found and even an enhance- 
ment is possible. This latter finding leaves the possibility 
open of substantial phonemic residual attractions between 
the quasiparticles competing with the residual repulsions 
in driving the system unstable toward long-wavelength 
charge instabilities. These two dynamical regimes are 
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also visible in the DMFT approach where the scale vf£[ 
is reflected in the width of the quasiparticle resonance. 
For frequencies smaller than this latter scale phonemic 
features are absent, while they are clearly present at high- 
energy in the Hubbard side-bands. 

The nearly static case co/vpq <C 1 again displays an 
intrinsic richness as far as momentum dependence is con- 
cerned: the strong suppression of g already occurring at 
small momenta becomes really very strong at large trans- 
ferred momenta. This suppression found both on general 
grounds within a Fermi-liquid analysis and within spe- 
cific field-theoretic treatments of the Hubbard-Holstein 
model 27, 28, 32, H, 4(| can account for the impressive 



elusiveness of phononic features in transport experiments 
in cuprates. Indeed the fact that the resistivity in the 
metallic phase does not display any clear phonon-related 
feature is naturally explained by the strong suppression 
of the e-ph coupling in transport processes, where very 
low-energy and large transferred momenta are involved. 
Thus the strongly correlated nature of the cuprates is 
the key ingredient to solve the puzzles related to the di- 
chotomic behavior of these materials, which display clear 
phononic features in some cases and none in others. 

All the above findings are only partly peculiar of the 
specific HH model, where the electron density locally in- 
volves both the e-e and e-ph coupling: The analysis of 
other models like the SSH essentially produce the same 
results 4(J 71 1. A qualitative difference only occurs for 
those phonons which couple to degrees of freedom, which 
are not severely stiffened by the proximity to a Mott- 
Hubbard phase. In this regard we reported the impor- 
tant case of Jahn- Teller phonons in the Fullerenes. It 
would be interesting to search for similar phononic (or 
even non-phononic) degrees of freedom in the cuprates. 
In this regard, the buckling modes in some cuprates are 
interesting candidates, which are presently being investi- 
gated in this perspective 72 1. 
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